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Abstract. Interference is one of the most fundamental features which characterizes 
quantum systems. Here we provide an exhaustive analysis of the interfere dynamics 
associated with wave-packet superpositions from both the standard quantum- 
mechanical perspective and the Bohmian one. From this analysis, clear and insightful 
pictures of the physics involved in this kind of processes are obtained, which 
are of general validity (i.e., regardless of the type of wave packets considered) 
in the understanding of more complex cases where interference is crucial (e.g., 
scattering problems, slit diffraction, quantum control scenarios or, even, multipartite 
interactions). In particular, we show how problems involving wave-packet interference 
can be mapped onto problems of wave packets scattered off potential barriers. 
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1. Introduction 

Over the last 15 years or so, fields such as the quantum information theory [1,2], 
quantum computation [1,2] and quantum control [3] have undergone a fast development. 
Simultaneously, due to the relevant role that entanglement plays in all these fields, the 
idea of this phenomenon being the most distinctive feature of quantum mechanics has 
also grown in importance, although this already infers from earlier works [4]. However, 
despite this relevance, in our opinion quantum interference (together with diffraction) 
can still be considered of more fundamental importance within the quantum theory, 
its most striking manifestation being, most surely, the two-slit experiment, which "has 
in it the heart of quantum mechanics. In reality, it contains the only mystery [. . .]", 
quoting Feynman [5]. Quantum interference is the direct, observable consequence of the 
coherent superposition of (quantum) probability fields. This is precisely what is special 
about it: interferences in quantum mechanics are not associated with or produced by 
a sudden transfer of energy (in the way of a perturbation) along a material medium, 
as happens with classical waves. It is worth stressing that, apart from its importance 
at a conceptual and fundamental level, quantum interference is also involved in a very 
wide range of experimental situations and applications. SQUIDs or superconducting 
quantum interference devices [6], the coherent control of chemical reactions [3], atom and 
molecular interferometry [7] (in particular, with BECs [8-10], where different techniques 
to recombine the split beams are used [11-14]) or Talbot /Talbot-Lau interferometry with 
relatively heavy particles (e.g., Na atoms [15] and BECs [16]) are just some examples. 
Moreover, it is also important to notice the role played by interference when dealing 
with multipartite entangled systems as an indicator of the loss of coherence induced by 
the interaction between the different subsystems, for instance. 

The actual understanding and interpretation of interference phenomena arises from 
the, somehow, physical "reality" ascribed to the superposition principle. According to 
this (mathematical) principle, wave fields can be decomposed and then recombined again 
to produce and explain the interference patterns observed. However, in Nature, wave 
fields constitute a whole and, therefore, it is important to determine what this behavior 
means when we are dealing with quantum phenomena displaying interference. Here 
we analyze such consequences in the case of coherent superpositions of Gaussian wave 
packets by means of properties which can be derived from standard quantum mechanics 
(but that, to our knowledge, are not considered at all in the literature) and also from 
Bohmian mechanics [17,18]. In the latter case, we will show that there exists a sort of 
mapping which allows us to translate wave-packet interference problems onto problems 
of single wave packets scattered off by potential barriers. The interest in this kind of 
superpositions arises from the fact that, rather than simple academical examples whose 
interference dynamics is well known in the literature [19,20], such superpositions are 
experimentally realizable (and, indeed, used) in atom interferometry [11-14]. In Fig. [TJ 
for instance, we have sketched two types of experiments where coherent superpositions 
of Gaussian wave packets can be produced. Moreover, it is also important to highlight 
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Figure 1. Two situations which give rise to a description of the wave function as a 
superposition of two wave packets at a certain time, (a) After reaching a beam-splitter, 
an incoming wave packet (1) splits into two ones (2). These wave packets move apart 
until they reach two reflectant mirrors, which redirect their propagation (3). The wave 
function describing the evolution of these two wave packets is a superposition which 
will display interference when they "collide" after some time in future (this fourth 
step has not been represented), (b) A wave packet confined in a potential well (1) 
is split into two wave packets by creating a barrier inside the potential (2) and then 
pulling outwards the two new potential wells (each containing a new wave packet). If 
the potential is turned off and the wave packets pulled inwards, the (superposition) 
wave function that describes them will again display interference (again, as in the 
previous case, we have not represented this step), as can be seen in BEC interferometric 
experiments. 

the interest that the analysis presented here can have to better understand scattering 
processes, quantum control scenarios or many-particle interactions, where interference 
plays a fundamental role, as well as to design, develop, improve and/or implement new 
trajectory-based algorithms [21], where one of the main drawbacks is the presence of 
nodes in the wave function produced by interferences. 

The organization of this paper is as follows. To be self-contained, in section [2] we 
provide a brief overview of quantum interference and Bohmian mechanics, as well as 
some interesting properties which will be used and further discussed in next section. 
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In section [3] we present the analysis of wave-packet superposition dynamics and the 
main results derived from it. Thus, we begin with a brief overview on (single) wave- 
packet dynamics and its properties, and then we introduce our non-standard analysis of 
two wave-packet interference processes and the analogy of this problem with wave- 
packet scattering off barriers. Finally, in section HI we provide a summary of the 
main conclusions in this work as well as a discussion on their interest in more complex 
problems. 



2. Interference and quantum trajectories 

Since the Schrodinger equation is linear, it satisfies the superposition principle. Thus, 
given ipi and ip 2 satisfying separately this equation, its superposition, 

tf(r,*) =il> 1 (r,t)+il> 2 (T,t), (1) 

will also be a valid solution. Now, the wave amplitude \I/ is not an observable magnitude, 
but the probability density, p = |\I/| 2 , which provides a statistical description of the 
system [22]. Due to the connection between \l/ and p, it is clear that the superposition 
principle does not hold for the latter. Note that, after expressing (i = 1, 2) in polar 
form, i.e., ipi = p] e %Si ^ h , we find 

p = pi + p 2 + 2 y /p 1 p 2 cos <p, (2) 

where <p = (S 2 — S\)/h. The interference term in Eq. (J2]) connects coherently the 
probability densities (pi and p 2 ) related to each partial wave (ipi and ip 2 , respectively). 

Since p describes the statistical distribution associated with a particle ensemble 
regardless of whether there is an interaction potential connecting them or not, one can 
define an associate probability current density as 

J = — Re[#*p^l = -— - , (3) 

m 2m 

which indicates the flow of such an ensemble, with p = —iftV being the momentum 
(vector) operator. Substituting ([I]) into (j3J), with the ipi expressed in polar form, yields 

j = i 

m 

which clearly shows that the superposition principle does not hold either for J. The two 
magnitudes p and J are related through the continuity equation, 

| + VJ = 0, (5) 

which can be easily derived from the Schrodinger equation after multiplying both sides 
by adding to the resulting equation its complex conjugate, and then rearranging 
terms. 

Within this hydrodynamical picture of quantum mechanics [23], one can always 
further proceed as in classical mechanics and determine the probability streamlines, i.e., 
the lines along which the probability flows. Consider the same polar ansatz used above, 



px VSi + p 2 VS 2 + v/pT^V^i + S 2 ) cosf + h {pT^pT - pT^pT) sin <p , (4) 
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but for the total wave function, \1/ = p 1 / 2 ^/^ After substituting it into the Schrodinger 
equation and rearranging terms, one obtains 

^+-V(pVS) = (7) 
at m 

from the real and imaginary parts of the resulting expression, respectively. Equation 
(I7j) is the continuity equation ([51), with J = (pVS)/m. Equation ([6]), more interesting 
from a dynamical viewpoint, is the quantum Hamilton- Jacobi equation, which allows us 
to reinterpret the whole quantum-mechanical formalism in terms of the eventual paths 
that the system can pursue when it is considered as a particle. These paths are defined 
as solutions of the equation of motion 

vs j ih 

j. _ _ (Q\ 

m p 2m 

where S represents a quantum generalized action. Note that considering Eq. ([8]) 
means introducing a new conceptual element into quantum mechanics: well-defined 
trajectories in both space and time. This new element constitutes the essence of what 
is known nowadays as Bohmian mechanics [17,18], a quantum mechanics based on the 
hydrodynamical Eqs. §6§ and ([7]) plus the particle equation of motion (JH]). Because of the 
relationship between particles and waves in Bohmian mechanics, the initial momentum 
of particles is predetermined by \l/(r, 0) and therefore it is not necessary to inquire about 
its value. Only initial positions are freely, randomly chosen, with the constraint that 
their distribution is given by p(r, 0). 

As infers from Eq. ([6]), the difference with respect to classical mechanics is that 
quantum trajectories evolve under the action of both the external potential V and the 
so-called quantum potential, 

h 2 V 2 p 1/2 

which introduces into the quantum motion the context-dependence and nonlocality 
necessary for the particles (distributed initially according to p(r, 0)) to reproduce 
the patterns of standard quantum mechanics when they are considered in a large 
statistical number. The quantum dynamics is thus ruled by a total effective potential 
V e ff(r, t) = V(r) + Q(r,t). An interesting property of the quantum potential can be 
found when one computes the energy expected value, 

E=(H) = (f) + (V), (10) 

which, in Bohmian mechanics, just consists of determining the ensemble average energy. 
The first term in (TTUT) is the expected value of the kinetic energy, with T = — (fr 2 /2m)V 2 
denoting the kinetic operator. When the polar ansatz of the wave function is considered, 
we obtain 

< f > = [ v2r - * + LI v( " vs)A -- < u > 
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From ([7]), we have 



VJdr = - I ^dr = -^- I pdr = 0, (12) 



dt dt 

since the probability density conserves in the whole space. In other words, if the system 
is closed, no probability can flow to or from the region where this system is defined; this 
is the probabilistic analog of the energy conservation principle. Therefore, (iTTj) and (fTUl) 
become 

and 

E = E k + Q + V = E k + Ves, (14) 

respectively, where E k is the average (expected value) contribution from the kinetic 
energy associated with each particle from the ensemble. Thus, notice that, although 
Q is usually assigned the role of a quantum potential energy, since it is obtained from 
the kinetic operator and contributes to the expected value of the kinetic energy, one 
can also consider it as a quantum kinetic energy. More importantly, putting aside such 
considerations, Q is the only responsible for making quantum dynamics so different from 
classical ones due to its nonlocal nature (E k and V only contain local information) — note 
that Q contains information about the whole quantum state and therefore supplies it to 
the particle at the particular point where it is located. From a practical or computational 
point of view, this explains why a good representation of the momentum operator has 
to contain a wide spectrum of momenta. Furthermore, (fT4"l) also indicates that classical- 
like behaviors appear whenever the part of the energy associated with the quantum 
potential becomes sufficiently small. 

To conclude this section, we are going to discuss another interesting property of 
quantum motion: in Bohmian mechanics two trajectories cannot cross the same point 
at the same time in the configuration space (note that, in classical mechanics, this only 
happens in the phase-space). In order to prove this, consider that two velocities are 
assigned to the same point r, v i(r) and t>2(r), with v\ 7^ v<i- From (jHJ) we know that these 
velocities are directly related to some wave functions ipi{ r ) an d V ; 2( r ), respectively. If 
both wave functions are solutions of the same Schrodinger equation, the only possibility 
for them to be different at the same space point x is that their phases differ, as much, 
in a time-dependent function (and/or space-independent constant), i.e., 

S 2 (T,t) = Sx(r,t)+(f>(t). (15) 

Now, since the velocities are the gradient of these functions, it is clear that v-i has to be 
equal to v\ necessarily if the wave function does not vanish at r. 
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3. Simple examples of quantum interference with Gaussian wave packets 

3.1. The free Gaussian wave packet 

For the sake of simplicity, we are going to consider one-dimensional Gaussian wave 
packets, although the results presented in this work can be generalized to more 
dimensions and other types of wave packets straightforwardly. The evolution of a free 
Gaussian wave packet can be described by 

where A t = (27raf )~ 1//4 and the complex time-dependent spreading is 
From (ITTl) . the spreading of this wave packet at time t is 




o~t = \<Tt\ = 

Due to the free motion, x t = x + v p t {y p = p/m is the propagation velocity) and 
E = p 2 /2m, i.e., the centroid of the wave packet moves along a classical rectilinear 
trajectory. This does not mean, however, that the expected (or average) value of the 
wave-packet energy is equal to E, which, according to Eq. (THl) . results 

E =^ + ^—. (19) 
2m Smal V ; 

In f )19p we thus observe two contributions. The first one is associated with the 
propagation or translation of the wave packet, while the latter is related to its spreading 
and has, therefore, a purely quantum-mechanical origin. From now on, we will denote 
them as E v and E s , respectively, noting that none of them depends on time. Expressing 
E s in terms of an effective momentum p s (i.e., E s = p 2 /2m), we find 

Ps = (20) 
2a 

which resembles Heisenberg's uncertainty relation. Note that, in the case of non- 
dispersive Airy wave packets [24], we will find p s = 0. 

With the previous definition, the wave packet dynamics can be referred to two well- 
defined velocities: v p and v s , whose relationship is shown to play an important role in 
the type of effects that one can observe when dealing with wave packet superpositions. 
To better understand this statement, consider Eq. f|T8|) . Defining the timescale 

r = 2mo-%/h, (21) 

we find that, if t <S t, the width of the wave packet remains basically constant with 
time, a t ~ o"o (i- e -> f° r practical purposes, it is roughly time-independent up to time t). 
On the contrary, if t ^> r, the width of the wave packet increases linearly with time 
(a t w ht/2mo- Q ). Of course, in between there is a smooth transition; from f|T8l) it is 
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shown that the progressive increase of a t describes a hyperbola when this magnitude is 
plotted vs time. These are "dynamical" relationships in the sense that the dynamical 
behavior of the wave packet is known comparing the actual time t with the effective 
time t. When we are dealing with general wave packets, a priori we do not have an 
effective time r to compare with. This inconvenience can be avoided if we consider the 
interpretive framework provided by v p and v s , which can be estimated from the initial 
state through Eq. ffl3|) . Within this framework, we can find a way to know which process 
— spreading or translational motion — is going to determine the wave packet dynamics, 
as follows. Expression (Tl~8l) can be rewritten in terms of v s as 



V V°o / 

Now, consider that t is the time that takes the wave-packet centroid to cover a distance 
d = v p t « a . Introducing this result into ff22l) . we obtain 



This time- independent expression shows that, only using information about the initial 
preparation of the wave packet, we can obtain information on its subsequent dynamical 
behavior. Thus, f < r is equivalent to having an initial wave packet prepared with 
v s <C v p : the translational motion will be much faster than the spreading of the wave 
packet. On the other hand, t 3> r is equivalent to v s ^> v p : the wave packet spreads very 
rapidly in comparison with its advance along x. A nice illustration of the "competition" 
between v p and v s can be found in [25], where (for |x | ~ °~o) it is shown that: (1) if v p 
is dominant, the asymptotic motion is basically a classical like motion (no spreading in 
comparison with the distances covered by particles), while (2) if v s is dominant, although 
the motion is classical-like (i.e., there is an asymptotic constant velocity, v s ), it has a 
purely quantum-mechanical origin and therefore eventual effects produced within this 
regime will also be purely quantum-mechanical. 

3.2. Dynamics of coherent wave-packet superpositions 

Depending on whether the interference is localized in time within a certain space region 
or it remains stationary after some time, we can distinguish two types of processes 
or experiments: collision-like (e.g., BEC interferometry) and diffraction-like (e.g., the 
double-slit experiment), respectively. These behaviors are associated with the ratio 
between v p and v s : collision-like experiments are characterized by v p ^> v s , while 
diffraction-like experiments by w p < f s . Now, consider a general coherent superposition 



where the partial waves ipi are normalized Gaussian wave packets, as in ffl6|) . which 
propagate with opposite velocities v p and where the centers are initially far enough in 
order to minimize the overlapping, i.e., pi(r, 0)p 2 (r, 0) ~ 0. As shown below, both the 




(22) 




(23) 



^ = Ci^l + C 2 ^ 2 , 



(24) 
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v p /v s ratio and the weighting factors c\ and c 2 are going to influence the topology of 
the quantum trajectories. Since the wave packets are well apart initially, we find that 
the weighting factors satisfy the relation |ci| 2 + | C2 1 2 = 1, which allows us to reexpress 
(I24j) more conveniently as 

* = c a (^i + v^j) , (25) 

where a = (c 2 /ci) 2 . From (1251) . we find that the probability density and the quantum 
current density read as 

p = c\ [px + ap 2 + 2y/a^/p 1 p 2 cos p] , (26) 
p{VS 1 + ap 2 VS 2 + Vay/pip 2 V(Si + S 2 ) cosip 



r 2 
m 



+ KyfaW'Vtf' - pT^pT) sm^j , (27) 

respectively. Also, dividing Eq. fl27|) by (1261) . as in (jSJ), we can obtain the associate 
quantum trajectories from the equation of motion 

1 piVSi + ap 2 VS 2 + y/a y /p 1 p 2 S/(Si + S 2 ) cosip 



m 



Pi + ap 2 + 2y/ay/pip 2 cosip 



k , 1/2^-, 1/2 1/2,-, 1/2 \ • 

+ v -^( Pl Vp 2 -p 2 V Pl )smy 
m pi + ap 2 + 2 v /a v /pip 2 cos^ 

As can be noticed, in this expression there are two well-defined contributions, which are 
related to the effects caused by the interchange of the wave packets (or the associate 
partial fluxes) on the particle motion (specifically, on the topology displayed by the 
corresponding quantum trajectories) after the interference process; note that the first 
contribution is even after interchanging only the modulus or only the phase of the wave 
packets, while the second one changes its sign with these operations. From the terms 
that appear in each contribution, it is apparent that the first contribution is associated 
with the evolution of each individual flux as well as with their combination. Thus, it 
provides information about both the asymptotic behavior of the quantum trajectories 
and also about the interference process (whenever the condition pi(r,t)p 2 (r,t) w is 
not satisfied). On the other hand, the second contribution describes interference effects 
connected with the asymmetries or differences of the wave packets. For instance, their 
contribution is going to vanish if they are identical and coincide at x = although their 
overlapping is nonzero. 

Let us first consider the collision-like case. For the sake of simplicity and without 
loss of generality, we will assume a — 1, with both wave packets being identical and 
propagating in opposite directions at the same speed. Also, for simplicity in the 
discussion, we will refer to the regions where ip\ and ip 2 are initially placed as I and 
II, respectively. In standard quantum mechanics, a physical reality is assigned to the 
superposition principle, which in this case means that, after the wave packets have 
maximally interfered at t™* x (panel (c) in the left column of Fig. [2]), ip\ moves to region 
II and ip 2 to region I. As mentioned above, the wave packets (or, more specifically, 
their associate probability densities) represent the statistical behavior of a swarm of 
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Figure 2. Left: From (a) to (d), snapshots illustrating the time-evolution of two 
colliding Gaussian wave packets, with v p = 10 and v s = 1. Center: The same process 
as in the left column, by visualized in terms of quantum trajectories. The perpendicular 
dashed lines mark cuts in time corresponding with the same labels as in the figures of 
the left column. Right: Interpretation of the process in the left column according to 
the quantum trajectories shown in the central column. 



identical, noninteracting quantum particles distributed accordingly. Therefore, within 
the standard view arising from the superposition principle, one would expect to observe 
crossings between trajectories in a certain time range around t^ax- However, this is not 
the behavior displayed by the Bohmian trajectories displayed in the central panel of 
Fig. [2j which avoid such crossings during the time range where interference effects are 
important. Note that quantum- mechanical statistics are characterized by keeping the 
coherence and transmitting it to the corresponding quantum dynamics, as infers from 
Eq. ( f28i) . Thus, the interference process has to be interpreted in a different way (than 
the standard quantum-mechanical one) when it is analyzed from a quantum trajectory 
perspective. As inferred from Eq. (f28i) . for identical (but counter-propagating) wave 
packets, the velocity field vanishes along x = at any time. This means that there 
cannot be any probability density (or particle) flux between regions I and II at any 
time. Therefore, trajectories starting in one of these regions will never cross to the 
other one and vice versa. Hence the final outgoing wave packets in panel (e) in the left 
column of Fig. [2] necessarily describe exactly the same swarms of particles associated 
with the wave packets that we had in their respective regions initially (see panel (a) 
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Figure 3. Bohmian trajectories associated with a Gaussian wave-packet superposition 
(grey) and a single Gaussian wave packet (red). As shown, the presence of the 
other wave packet makes the dynamics to avoid the crossing in the central part; the 
asymptotic part remains being the same as without presence of other wave packet. 
Here v„ — 10 and v s — 1. 



of the same figure). The whole process can be then understood as a sort of bouncing 
motion of the wave packets once they have reached the intermediate position between 
them, as schematically represented in the right column of Fig. El 

But, if the two swarms of trajectories do not cross each other, why the swarm 
associated initially with one of the wave packets behaves asymptotically as associated 
with the other one, as seen in Fig. [3P According to the standard description based on the 
superposition principle, the wave packets cross. This can be understood as a transfer or 
interchange of probabilities from region I to region II and vice versa. On the other hand, 
from a dynamical (quantum trajectory) viewpoint, this also means that the sign of the 
associate velocity field will change after the collision. That is, before the collision its 
sign points onwards (towards x = 0) and after the collision it points outwards (diverging 
from x = 0) — at t™^ x it does not point anywhere, but remains steady. Thus, like in 
a particle-particle elastic scattering process particles exchange their momenta, here the 
swarms of particles will exchange their probability distributions "elastically" . This is 
nicely illustrated in Fig. [3j after the collision, the two swarms of trajectories (grey lines) 
bounce backwards and follow the paths that would be pursued by non deflected particles 
(red lines). We can describe this process analytically as follows. Initially, depending on 
the region where the trajectories are launched from, they are described approximately 



by 



ri ~ VSi/m 



or 



r n ~ VS 2 /m. 



(29) 



Now, asymptotically (for t 3> £™a X ), Eq. (1281) can be expressed as 

1 P1VS1 + P2VS2 
r « ■ , 

m pi + p 2 



(30) 
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where we make use of the approximation pi(r, t)p 2 (r, t) ~ 0. Note that this 
approximation also means that the total probability density, p = p\ + p 2 , is nonzero 
only on the space regions covered by either p\ or p 2 . Specifically, in region I we will 
have p ~ P2 and in region II, p « p 1 . When this result is introduced into Eq. (130]) . we 
finally obtain 

ii ~ VS 2 /m and rn ~ VS'i/m, (31) 

which reproduce the dynamics observed asymptotically. This trajectory picture thus 
provides us with a totally different interpretation of wave packet interference with respect 
to the standard one: although probability distributions transfer, particles remain always 
within the domains defined by their corresponding initial distributions (which have been 
here as regions I and II). The non- crossing property of Bohmian mechanics (see section[2]) 
thus makes apparent the constraint that exists for the quantum probability flux, which 
goes beyond the separability of fluxes implicit in the superposition principle. 

In the diffraction-like case, the spreading is faster than the propagation, this being 
the reason why we observe the well-known quantum trajectories of a typical two-slit 
experiment in the Fraunhofer region [19,26]. It is relatively simple to show [20] that, in 
this case, the asymptotic solutions of Eq. ( 1281) are 

x(t) « 2vm — (v s t) , (32) 
x 

with n = 0, ±1,±2, ... That is, there are bunches of quantum trajectories whose 
slopes are quantized quantities (through n) proportional to v s . This means that, when 
Eq. (|2"5|) is integrated exactly, one will observe quantized bunches of trajectories which, 
on average, are distributed around the value given by Eq. fl32l . 

3.3. Asymmetric coherent superpositions 

In the previous subsection we have studied interference processes with identical Gaussian 
wave packets propagating symmetrically with respect to x — 0. However, what happens 
in other more general situations where neither the wave packets are identical nor their 
weights? Although one could think of many different possibilities, it is enough to 
consider three cases in order to already obtain an insight on the related physics. The 
criterion followed to classify these cases is based on which property of the wave packets 
is varied (with respect to the symmetric case described in section 13. 2j) or considered at 
a time: 

• Case A: Different modulus of the initial average momentum (|poi| 7^ bed)- 

• Case B: Different initial spreading (o"oi ^ o~q 2 ) . 

• Case C: Different weights (ci 7^ c 2 or, equivalently, a^l). 

Since the collision-like case is simpler to understand than the diffraction-like one, below 
we will assume v p 3> v s , although generalizing to any v p /v s ratio is straightforward. 
Moreover, as before, the initial distance between both wave packets is such that the 
initial overlapping is negligible. 
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Figure 4. (a) Probability density at t — (black-dashed line) and t — 0.6 (blue-solid 
line) for a wave packet superposition with p i — 10 and p 2 = —30. The value of the 
other relevant parameters are <7oi = C02 = 0.5 and a = 1. (b) Quantum trajectories 
associated with the dynamics described by the situation displayed in panel (a). 



Case A is represented in Fig. HI where we can observe that varying the moduli of 
the initial average momenta only produces an asymmetric shift of the final positions 
of the wave packets with respect to x = (see panel (a)), which is a result of the 
distortion of the boundary between regions I and II (the final wave packets are not 
symmetrically distributed with respect to x = 0). This boundary evolves in time at the 
same constant velocity than the corresponding expected value of the velocity for the 
total superposition, 

Pa (*1IpI*a) P01+P02 /qq\ 

V A = = = 7, > 

because both wave packets are identical and normalized (see Eq. ( fl6l) ). In particular, 
since poi = 10 and P02 = —30, we have va = —10. Integrating the equation of motion 
xa = va (according to Ehrenfest's theorem [27]), we obtain 

x A {t) = ^0 + v A t, (34) 

with xq = (ty*\x\ty) = (x i +Xo 2 )/2. Equation (134)) describes the time-evolution of the 
boundary between regions I and II in this case — note that for p 02 = —poi, we recover 
again the time-independent boundary found for identical wave packets, studied in the 
previous section. As expected, this boundary also defines the non-crossing line for the 
associate quantum trajectories, which forbids the trajectory transfer from one region 
to the other one and vice versa, as seen in Fig. H^b). Because of this property, we 
find that the two wave packets behave like two classical particles undergoing an elastic 
scattering: there is only transfer (indeed, exchange) of momentum during the scattering 
process, but the net balance of probability is zero since no trajectories are transferred. 
Moreover, note that indeed this effect will not be noticeable unless one looks at the 
quantum trajectories — as seen in panel (a), the evolution of the wave packets does not 
provide any clue on it. 
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Figure 5. (a) Probability density at t — (black-dashed line) and t — 0.8 (blue-solid 
line) for a wave packet superposition with cr 01 = 0.5 and 002 = 1-5. The value of the 
other relevant parameters are \poi\ = \po2\ = 20 and a = 1. (b) Quantum trajectories 
associated with the dynamics described by the situation displayed in panel (a). 



In case B, since both velocities are equal in modulus, after interference the centers 
of the outgoing wave packets occupy symmetric positions with respect to x = 0, as 
seen in Fig. |Sfa) (only the different width disturbs the full symmetry). However, as in 
case A, the lack of total symmetry also causes the distortion of the boundary between 
regions I and II (see panel (b)), with a similar behavior regarding the non-crossing 
(or non-transfer) property. Now, since the modulus of both momenta are the same, 
(^ r *|p|^ r ) = and therefore we cannot appeal to the same argument as before to explain 
this distortion effect. However, there is still a sort of effective "internal" momentum 
which depends on the ration between the widths of the wave packets. Assuming that 
the spreading is basically constant for the time we have propagated the trajectories, the 
corresponding effective velocity is given by 

vb = _! -i = - — Vq, (35) 

with Voi = vq = —vq2 {vq > 0); since the width of the wave packets varies in time, a 
slight time- dependence can be expected in vb, although in our case, as can be noticed 
from Fig. [5](b), it can be neglected for practical purposes. Substituting the numerical 
values used in the propagation into ( |35i) . we find vb = 5. The boundary or non-crossing 
line will be then given by 

^b(*) = + v B t, (36) 

where 

gpl 1 Xqi + CTpa 1 XQ2 f 0-Q2-0-Ql \ , , 

xo- =5— ■ — — = -[- — — — x o, [o<) 

a 0l + a 02 \ a 01 + &02/ 

with X02 = —xoi = xo (xo > 0) — note that, otherwise, xo = 0, as infers from Fig. [5(a), 
since the larger width of one of the wave packets balances the effect of the larger height 
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of the other one. Despite this internal redistribution or balance of momentum, it is clear 
that we find again, as in case A, an elastic collision-like behavior. 

The effect of an effective internal momentum could be explained by considering 
that, in this case, the relative size of the wave packets acts like a sort of quantum inertia 
or effective mass. Consider the trajectories associated with the wave packet with smaller 
value of cr (black trajectories in Fig. E(b)). As can be noticed, it takes approximately 

0. 2 time units the whole swarm to leave the scattering or collision region (along the 
non-crossing boundary) — or, equivalently, to revert the sign of the momentum of all 
the trajectories constituting the swarm and get their final asymptotic momenta. On 
the other hand, the trajectories associated with the wave packet with larger <7o (red 
trajectories in Fig. E(b)) revert their momenta much faster, in about 0.1 time units 
— of course, we are not considering here the time that trajectories keep moving along 
the boundary, since the total "interaction" time has to be the same for both swarms 
of trajectories. Thus, we find that the larger the spreading momentum, the larger the 
quantum inertia of the swarm of particles to change the propagation momentum and, 
therefore, to reach the final state. 

We would also like to note another interesting property associated with case 
B. Consider that both wave packets lack their corresponding normalizing factor A t 
when introduced in the superposition. If we then compute the expected value of the 
momentum, for instance, we obtain 

(*£|p|* B ) 0-01^01 + 0-02^02 fa 01 -a 02 \ , . 

iP)B ~ " o 01 + a 02 ~ V^T + ^J P °' (38) 

1. e., there should be a certain "drift" towards region I, such as in case A — and the 
same holds if we compute instead the expected value of the position. Note that in 
the previous case the normalization of each Gaussian wave packet produces a balance: 
the probability with which each wave packet contributes to the superposition is the 
same (c\ = c\ = 1/2) because the width of one compensates the height of the other, as 
explained above. Therefore, the expected value of both position and momentum have to 
be zero. However, this compensation does not happen now: both wave packets have the 
same height although their widths differ, thus contributing with different probabilities 
P to the superposition, 

P 1 = am and P 2 = ° m , (39) 

O i+O 2 O 01 +cr 02 

which produce the results observed in Fig. [6] (again, we assume that a t ~ o"o for the 
time considered). However, by inspecting (J3S|) . we note that if we add the averaged 
momentum p-Q = m%, the total average momentum vanishes. Somehow the averaging 
defined by (j38p acts as in classical mechanics, when a certain magnitude (e.g., the 
position or the momentum) is computed with respect to the center of mass of a system. 
Here, ps is the magnitude necessary to reset the superposition of non normalized 
Gaussian wave packets to a certain "center of spreading" . On the other hand, it is also 
important to stress the fact that, in this case, the clear boundary between the swarms 
of trajectories associated with each initial wave packet disappears. Now, although there 
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Figure 6. Bohmian trajectories associated with a Gaussian wave-packet superposition 
with ctoi = 0.5 and 002 = 1.5. The value of the other relevant parameters are 
P01 = — P02 = 20 and a = 1. 



is still a boundary, it does not prevent for the transfer of trajectories from one region 
to the other, as before. This effect, similar to consider inelastic scattering in classical 
mechanics, arises as a consequence of having wave packets with different probabilities, 
which we analyze below. 

As seen above, unless there is an asymmetry in the probabilities carried by each 
partial wave in the superposition, there is always a well-defined boundary or non-crossing 
line between regions I and II. However, in the second example of case B, we have observed 
that it is enough an asymmetry in the probability distribution to break immediately 
the non-crossing line. Instead of using non normalized Gaussian wave packets, let 
us consider case C, which is equivalent although the asymmetry is caused by a / 1 
instead of the wave-packet normalization. As we have seen in the two previous cases, 
provided both wave packets have the same weight (or, at least, both contribute equally 
to the superposition), the corresponding quantum trajectories, even if they are evenly 
distributed (equidistant) along a certain distance initially, they are going to give a 
good account of the whole dynamics. However, if the weights change, the same does 
not hold anymore. Somehow equal weights (or equal probabilities) means the division 
of the coordinate space into two identical regions, each one influenced only by the 
corresponding wave packet. In other words, there are always two well-defined swarms 
of trajectories, each one associated with one of the initial wave packets. When the 
weights change, apparently we have something similar to what we have been observing 
until now: the wave packets exchange their positions (see Fig. [3(a)). However, when 
we look at the associate quantum trajectories, we realize that there is transfer or flow 
of trajectories from one of the initial swarms to the other one. This transfer takes 
place from the swarm with larger weight to the lesser one, thus distorting importantly 
the boundary between regions I and II, as seen in Fig. [3(b), where this boundary lies 
somewhere between the two trajectories represented in blue (the trajectories that are 
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Figure 7. (a) Probability density at t = (black-dashed line) and t = 0.8 (red-solid 
line) for a wave packet superposition with a = 0.5 (i.e., c\ = 2/3 and c| = 1/3). The 
value of the other relevant parameters are poi = —P02 = 10 and aoi = 002 = 0.5. 
(b) Quantum trajectories associated with the dynamics described by the situation 
displayed in panel (a). The initial conditions for the quantum trajectories have been 
assigned by considering the different weights associated with each wave packet. 



closer to x = in each swarm). However, note that this does not imply that the 
number of trajectories varies in each region asymptotically, but only the number of 
them belonging to one or the other wave packet. Thus, if initially we have Ni oc c\ 
trajectories associated with ipi and N 2 oc c| with ip 2 , asymptotically we will observe 
N[ = (l — a)Ni oc c 2 and N' 2 = N 2 + aN 1 oc c\ due to the trajectory transfer. Unlike the 
two previous cases discussed above, this process can be then compared with inelastic 
scattering, where, after collision, not only the probability fluxes but also the number 
of particles changes. Accordingly, it is also important to mention that, due to the 
trajectory transfer, representations with evenly distributed trajectories are not going 
to provide a good picture of the problem dynamics, as infers from Fig. [61 Rather, we 
can proceed in two different ways. The obvious procedure is to consider initial positions 
distributed according to the corresponding (initial) probability densities. This procedure 
carries a difficulty: the number of trajectories needed to have a good representation 
of the dynamics may increase enormously depending on the relative weights. The 
second procedure is to consider evenly spaced values of the probability density and 
allocate at such space points the initial positions of the trajectories, as we have done to 
construct Fig. [3(b). In this case, although the trajectories will not accumulate exactly 
along the regions with larger values of the probability density, this construction has 
the advantage that we can follow the transport of equi-spaced probabilities along each 
particular trajectory. To make more apparent the difference in the relative number 
of trajectories (probability density) associated with each initial wave packets, we have 
considered N x = 23 and N 2 = 11 (i.e., N 2 /N 1 ~ 0.48 w c\jc\ = 0.5). After scattering, 
the numbers that we have are N[ = 11 and N' 2 = 23, which are in the expected ratio 
c 2I c 1 = C l/ C 2 = 2- 
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As said at the beginning of this subsection, here we have only considered the 
collision-like case. The same kind of results are expected for the analogous diffraction- 
like cases, with the addition that they manifest as a loss of fringe visibility. Thus, for 
cases A and B, one could appreciate a well define non-crossing line and a loss of fringe 
visibility due to the divergent velocities, v p and v s , respectively. And, for case C (or the 
second example in case B) , the loss of fringe visibility would be caused by the transfer 
of trajectories, which would also lead to the distortion of the non-crossing line. 

3.4- Coherent superpositions, potential barriers and resonances 

Finally, here we are going to analyze an interesting effect associated with the non- 
crossing property described in the previous sections. Consider a Gaussian wave packet 
scattered off an impenetrable potential wall (for simplicity, we will assume v p > v s for 
now). After some time the wave packet will collide with the wall and then part of it will 
bounce backwards. The interference of the forward (/) and backward (b) wave packets 
will lead to a fringe-like pattern similar to those observed in the previous sections (see 
red-solid line in Fig. 0(a)), with also a time, t™* x , for which the interference fringes are 
maximally resolved. Putting aside the initial Gaussian shape of the wave packet (and, 
therefore, the effects associated with v s ), if this process is represented as 

m = ^ f + ^ b ~ e imv v x/h + e -imv p x/h (40) 

at t^ax) we obtain p(x) ~ cos 2 (mv p x/h). The distance between two consecutive minima 
is then w Q = irh/mv p , which turns out to be the same distance between two consecutive 
minima in the two wave-packet interference process (see black-dashed line in Fig. El(a)). 
That is, although each process has a different physical origin (barrier scattering vs wave 
packet collision), the effect is similar — there is a certain shift in the position of the 
corresponding maxima (~ vr/2), which arises from the fact that, in the case of barrier 
scattering, the impenetrable wall forces the wave function to have a node at x = 0. If 
now we go to the corresponding quantum trajectories, we observe (see Fig. [51(b), with 
red-solid line) that as the wave packet starts to "feel" the presence of the wall, the 
trajectories bend gradually (in the x vs t representation) and then start to move in the 
opposite direction. When these trajectories are compared with those associated with 
the problem of the two wave-packet superposition, the resemblance between trajectories 
with the same initial positions is excellent, except in the interference region due to the 
different location of the nodes of the corresponding wave functions — these differences 
are the trajectory counterpart of the shift mentioned before. 

From the previous description one might infer that, since the dynamics for x < 
and for x > do not mix (due to non-crossing), each half of the central interference 
maximum arises from different groups of trajectories. Thus, in principle, one should be 
able to arrange the impenetrable wall problem in such a way that allows us to explain 
this effect. Within this context, although all the peaks have the same width as in the 
wave-packet superposition problem, the closest one to the wall should have half such a 
width, i.e., w ~ nh~/2p = wq/2. Due to boundary conditions and the forward-backward 
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Figure 8. (a) Probability density at t = 0.3 for the collision of a Gaussian wave packet 
off: an external impenetrable wall (red-solid line) and another (identical) Gaussian 
wave packet (black-dashed line). To compare, the maxima of both probability densities 
are normalized to unity, (b) Bohmian trajectories illustrating the dynamics associated 
with the two cases displayed in part (a). Here, v p — 10 and v s = 1. 



interference discussed above, it is clear that this peak cannot arise from interference, 
but from another mechanism: a resonance process. Therefore, apart from the wall, we 
also need to consider the presence of a potential well. In order to observe a resonance or 
quasi-bound state, the width of this well should be, at least, of the order of the width 
w of the bound state. From standard quantum mechanics, we know that in problems 
related to bound states in finite well potentials the relationship [27] 

V a 2 = n^- (41) 

always appears, where a is the half- width of the well (a = w/2) and n is an integer 
number. The eventual solutions (bound states) are then observable or not depending 
on whether the condition which it might correspond will be in consonance or not with 
this condition. In our case, we can use ( ]4~T1) to obtain an estimate of the well depth, 
which results 

16 p 2 

V o = ~2&- ( 42 ) 
ir z 2m 

when we assume n = 1. Now, we have then a potential which presents a short-range 
attractive well before reaching the impenetrable wall, 

x < —w 

V(x) = { -V -w < x < . (43) 

oo < x 

If we compute now p(r,t) at t^ x , we obtain the result displayed in Fig. [9](a). As can 
be noticed, now there is an excellent matching of the peak widths, with the closest one 
to the wall being half-width when compared with the remaining ones — the associate 
quantum trajectories are displayed and compared in Fig. Mjo). 
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Figure 9. (a) Probability density at t = 0.3 for the collision of a Gaussian wave packet 
off: the external potential described by Eq. (|43p (red-solid line) and another (identical) 
Gaussian wave packet (black-dashed line). To compare, the maxima at x ~ —0.32 of 
both probability densities are normalized to the same height (the maximum at x — of 
the two wave packet probability density being set to unity), (b) Bohmian trajectories 
illustrating the dynamics associated with the two cases displayed in part (a). Inset: 
enlargement of the plot around t = 0.15 to show the action of the potential well on the 
trajectories started closer to x — 0. Here, v p — 10 and v s = 1. 



The equivalence between the two wave packet collisions and the scattering of a 
wave packet off a potential is not restricted to the condition v p > v s , but it is of general 
validity. As shown in Fig.[5j it also holds for the diffraction-like situation, i.e., v p < v s . In 
Fig. ITDT a) . we show that half of the diffraction-like pattern is again well reproduced after 
replacing one of the wave packets by an external potential, and the same also happens for 
the corresponding quantum trajectories (see Fig. fTOT b)). However, in order to find these 
results, now a subtlety has to be considered: the central diffraction maximum increases 
its width with time. In terms of simulating this effect with a potential function, it is 
clear that the width of the potential well should also increase with time. Thus, we need 
to consider a "dynamical" or time-dependent potential function rather than a static 
one, as done before. In order to determine this potential function, we proceed as before. 
First, we note that the two wave-packet collision problem, Eq. ( l24l) . is explicitly written 
in terms of Gaussian wave packets (see Eq. (fl6|) ) as 

\J/ e -{x+x t ) 2 / Aa t cro+ip{x+xt) / h+iEt/ h _|_ ^-(x-xt) 2 / ' Aa t m)-ip{x-x t ) / 'fl+iEt / 'h (44) 

where 

x t = x - v p t (45) 

(for simplicity, we have neglected the time-dependent prefactor, since it is not going 
to play any important role regarding either the probability density or the quantum 
trajectories). The probability density associated with (T44T) is 

P (X, t) ~ e -ix+*l)M + e -0»-s?)/2o? + 2e -(x 2 +x 2 )/2a 2 ^ y ^ ? ^ 
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Figure 10. (a) Probability density at t = 5 for the collision of a Gaussian wave 
packet off: the external, time-dependent potential described by Eq. (|53[) (red-solid line) 
and another (identical) Gaussian wave packet (black-dashed line). To compare, the 
maxima at x w —5 of both probability densities are normalized to the same height (the 
maximum at x — of the two wave packet probability density being set to unity) . (b) 
Bohmian trajectories illustrating the dynamics associated with the two cases displayed 
in part (a). Here, v p = 0.1 and v s = 1. 
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As can be noticed, (]4"6l) is maximum when the cosine is +1 (constructive interference) 
and minimum when it is —1 (destructive interference). The first minimum (with respect 
to x = 0) is then reached when f(t)x = it, i.e., 
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(49) 



which basically is zero if the initial distance between the two wave packets is relatively 
large when compared with their spreading. 

In Fig. HIT a), we can see the function x min (t) for different values of the propagation 
velocity v p . As seen, x m i n (t) decreases with time up to a certain value, and then increases 



again, reaching a linear asymptotic behavior. From 
of x min (t) is reached at 



4mo"o 



h 2 xo 




8l) we find that the minimum value 



(50) 



The linear time-dependence at long times is characteristic of the Fraunhofer regime, 
where the width of the interference peaks increases linearly with time. On the other 
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Figure 11. (a) Plot of x m - m as a function of time for different values of the propagation 
velocity: v p = 0.1 (black solid line), v p = 2 (red dotted line), v p = 10 (green dashed 
line) and v p = 100 (blued dash-dotted line), (b) Plot of Vq as a function of time for 
the same four values of v p considered in panel (a). In all cases, v s = 1. 



hand, the fact that, at t = 0, x m i n (t) increases as v p decreases (with respect to v s ) could 
be understood as a "measure" of the coherence between the two wave packets, i.e., how 
important the interference among them is when they are far apart (remember that, 
despite their initial distance, there is always an oscillating term in between due to their 
coherence [28]). Note that this is in accordance with the standard quantum- mechanical 
arguments that interference-like patterns are manifestations of the wavy nature of 
particles, while scattering-like ones display their corpuscle nature (more classical-like). 
Thus, as the particle becomes more "quantum-mechanically" , the initial reaching of 
the "effective" potential well should be larger. And, as the particle behaves in a more 
classical fashion, this reaching should decrease and be only relevant near the scattering 
or interaction region, around x = 0. From Eq. (|48|) . two limits are thus worth discussing. 
In the limit p ~ 0, 

x t 

and t min ~ r. In the long-time limit, this expression becomes x min (t) w (7ih/2m)(t/x ), 
i.e., x m ; n increases linearly with time, as mentioned above. On the other hand, in the 
limit of large o"o (or, equivalently, v p ^> v s ), 

x min (t) w ^ (52) 
2p 

and t m i n ~ 0. That is, the width of the "effective" potential barrier remains constant in 
time, this justifying our former hypothesis above, in the scattering-like process, when 
we considered w ~ nh/2p. 
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After Eq. f)48p . the time- dependent "effective" potential barrier is denned as 

{0 x < x min (t) 

-V [x min (t)] x min (t)<x<0 . (53) 

oo < x 

with the (time-dependent) well depth being 

2H 2 1 

V [x min (t)} = ^— . (54) 

m x 2 miQ {t) 

The variation of the well depth along time is plotted in Fig. [TlTa) for the different values 
of v p considered in Fig. [TTT b). As seen, the well depth increases with v p (in the same 
way that its width, x m { n , decreases with it) and decreases with time. For low values of 
v p , there is a maximum, which indicates the formation of the quasi-bound state that will 
give rise to the innermost interference peak (with half the width of the remaining peaks, 
as shown in Fig. [T0l(a)). Note that, despite the time-dependence of the well depth, in 
the limit v p ^> v s , we recover Eq. ()42l) . 

We have shown that the problem of the interference of two colliding wave packets 
can be substituted by the problem of the collision of a wave packet off a potential barrier. 
Although the model that we have presented is very simple, it is important to stress that 
it reproduces fairly well the dynamics involved — of course, one can always search for 
more refined and precise models. The fact that one can make this kind of substitutions 
puts quantum mechanics and the superposition principle on the same grounds as two 
important frameworks in classical physics (one could think of many other situations, but 
these two ones are particularly general and well known). The first one is widely used 
in classical mechanics (and then, in its quantized version, also in quantum mechanics): 
it is the equivalence between a two-body problem and a one-body problem acted by a 
central force. As can be noticed, our reduction is of the same kind: a two wave-packet 
dynamical process can be reduced to the dynamics of a single wave packet subjected to 
the action of an external potential. The second framework is the one provided by the so- 
called method of images, widely used in electrostatics: the interaction between a charge 
distribution and a conductor can be replaced by the interaction between such a charge 
distribution and another virtual one. This would be the reciprocal situation to ours: the 
dynamics induced by an external potential on a wave packet can be translated as the 
dynamics taking place when two wave packets (ours and a virtual one) are considered. 



4. Final discussion and conclusions 



Here we have presented a detailed analysis of quantum interference within the framework 
of quantum trajectories. As we have shown, there are at least two important properties 
which are masked when interference is studied from the standard point of view, but 
that become apparent as soon as we move to the Bohmian domain: (a) it is always 
possible to determine uniquely the departure region of a particle from the final outcome 
and (b) interference can be described as an effect similar to that of particle-particle 
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Figure 12. Schematic representation of the two types of processes that may take 
place in two wave-packet collisions: (a) elastic-like scattering and (b) inelastic-like 
scattering. See text for details. 

collisions. Bohmian trajectories allow one to determine the origin of the particles that 
contribute to each final wave packet, contrary to the standard view of interference 
processes, where the final wave packets are commonly associated with the time-evolved 
initial ones, according to the superposition principle. Nevertheless, as has also been 
shown here, putting aside the superposition principle and considering the properties of 
the quantum current density, it is also possible to reach similar conclusions in standard 
quantum mechanics (although, to our knowledge, one cannot find this in the literature). 
Of course, although we might know that quantum particles are associated with one 
wave packet or the other during the whole propagation, this should not be mistakenly 
considered as to have a complete knowledge of their corresponding initial positions in 
the experiment — remember that the initial positions distribute in a random fashion 
following the initial probability density. 

To some extent, the property (a) can then be regarded as a sort of distinguishability, 
which should not be confused with that one arising when one considers the problem 
of identical (fermion or boson) particles. Rather, this issue is connected with the 
property (b) in the sense that it allows us to describe the process, in the collision- 
like CclSC ; clS classical collisions. Two situations can then happen. If the probability is 
equi-distributed between the two wave packets, the effect is similar to classical elastic 
particle-particle collisions, where only the momentum is transferred. In our case, what 
is transferred is the probability distribution through the probability density current, 
since the number of particles associated with each wave packet remains the same. This 
situation, represented schematically in Fig. fl~2l(a). corresponds to the cases A and the 
first one of B described in section 13.31 On the other hand, when the probabilities 
are asymmetrically distributed, there is a probability transfer which translates into a 
trajectory transfer — the well defined non-crossing line then disappears. This behavior, 
which is similar to a classical inelastic particle-particle collision, is displayed in Fig. [T27b) 
and corresponds to the case C and the second example of B (see section l3~3l) . Rather than 
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being just a conceptual discussion, it is worth noticing that this analysis can be of much 
applied importance in the design, development, improvement and/or implementation of 
new trajectory-based algorithms [21], where one of the main problems is precisely the 
way to deal with separate (in space) wave packets and the nodes emerging from their 
interference. 

As has been noticed, observing collision-like behaviors or not depends on the 
relation between the spreading and the propagation velocities of the wave packets. If v p is 
larger than v s , one will observe collision-like behaviors, which, to some extent, is the same 
to say that the wave packets behave like classical particles or corpuscles. Of course, note 
that when we use the concept "collision" we are not referring exactly to a true particle- 
particle collision, since the two wave packets indeed refer to the same system (which, 
by means of any of the procedures mentioned in section [H for instance, has been split 
up and converted into a so-called Schrddinger cat system). On the other hand, when v s 
dominates the dynamics, the behavior is diffract ion- like. The wave packets then behave 
in a wavy manner and diffraction patterns appear after their interference. 

Finally, it is also very important the fact that wave-packet interference problems can 
be understood within the context of scattering off effective potential barriers. As is well 
known, in classical mechanics one can substitute a particle-particle scattering problem 
by that of an effective particle (associated with the center of the mass of the particle- 
particle system) acted by an (also effective) interaction central potential. Similarly, here 
we have shown that any interference process (either collision-like or diffraction-like) can 
also be rearranged in such a way that the two wave-packet problem can be replaced 
by that of a (single) wave-packet scattering off an effective (time-dependent) potential 
barrier. In this regard, we would like to note that this property could be considered a 
precursor of more refined and well-known methods used to deal with many body systems, 
where the many degrees of freedom are replaced by a sort of effective time-dependent 
potential [29]. Furthermore, it is worth stressing that, in order to describe properly 
the interference fringes during the process, we have shown that these potentials have to 
support temporary bound states or resonances. In those cases where the non-crossing 
boundary disappears because the probability is not equally distributed between the two 
wave packets, the impenetrable barriers should be replaced by "transparent" or step 
barriers. Thus, in the same way that after the collision process part of the trajectories 
associated with one of the wave packets become attached to the other wave packet, 
after the collision with the transparent barrier part of the trajectories will move on the 
barrier. 
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